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Abstract. Although often chosen because of simplicity, a single scalar field does not 
provide a general parametrization of an MSSM flat direction. We derive a formalism 
for a class of gauge invariant polynomials which result in a multifield description of 
the flat directions. In contrast to the single field case, the vanishing of the gauge 
currents yields an important dynamical constraint in the multifield framework. We 
consider in detail the example of the H^L flat direction and study the dynamical 
evolution during and after inflation. We highlight the differences between the single 
and the multifield flat directions. We show that in the multifield case the field space 
has an intrinsic curvature and hence unsuppressed non-minimal kinetic terms for the 
flat direction scalars arise. Also the phases of the individual components evolve non- 
trivially right after inflation, charging the components of the condensate and producing 
an enhanced entropy after the decay of the condensate, which is due to cross-coupling 
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of different lepton flavours in the F term. However, the quahtative features of the 
single field Affleck-Dine baryogenesis, such as the produced total charge, remain largely 
unchanged. 

1. Introduction 

Super symmetric gauge theories come with gauge invariant polynomials along which 
the scalar potential vanishes at a classical level PJ. These correspond to the flat 
directions in the moduli space of the scalar fields. Within the Minimal Supersymmetric 
Standard Model (MSSM) the classical potential is given by the sum of the F and D term 
contributions, which vanish individually along these fiat directions. Flatness is however 
lifted by soft supersymmetry breaking, but a non-renormalization theorem guarantees 
that the fiat direction does not obtain any renormalizable perturbative superpotential 
correction However, non-perturbative (supergravity- induced) corrections will also 
lift the flatness. 

MSSM fiat directions can give rise to a host of cosmologically interesting dynamics 
(for a review, see pj). These include Affleck-Dine baryogenesis jHIHlEl, the cosmological 
formation and fragmentation of the MSSM fiat direction condensate and subsequent 
Q-ball formation [3 IH1 El reheating the Universe with Q-ball evaporation P^ . 
generation of baryon isocurvature density perturbations fT^, as well as curvaton 
scenarios where MSSM flat directions reheat the Universe and generate adiabatic density 
perturbations [THj. Adiabatic density perturbations induced by fluctuating inflaton- 
MSSM flat direction coupling has also been discussed [Hj . 

In the MSSM all the flat directions have been parameterized by gauge invariant 
monomials [HIE]- In such a case the MSSM flat direction is some linear combination of 
the MSSM scalars and can be thought of as a trajectory in the moduli space described 
by a single scalar degree of freedom. A simple example is the HuL flat direction 0, 
defined as {Hu)i = L2 = (indices here label SU(2)). However, a flat direction can 
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also be constructed as a gauge invariant polynomial, where the previously mentioned 
monomials generate the polynomial. This means that the fiat direction is a non-linear 
combination of the MSSM fields. It no longer can be described by a single scalar field. 
Such a gauge invariant polynomial, rather than a monomial, is needed e.g. when for 
HuL instead of a single leptonic generation one considers simultaneously all three. 

The purpose of the present paper is to illustrate the dynamics of Affleck-Dine 
condensate described by a gauge invariant polynomial. Previously Senami and 
Yamamoto jTHj have considered such a case in MSSM and its extensions but they used 
an extrinsic point of view without gauge invariant polynomials. We will discuss the 
dynamics from an intrinsic point of view i.e. using coordinates that describe only the 
behaviour of the flat direction. Given that the fiat direction is the minimum of energy 
the two approaches are equivalent. 

We will also discuss how a non-renormalizable superpotential correction can lift such 
a multi-field fiat direction. We will highlight the formation and charge conservation of 
the condensate along with the gauge conditions required to be satisfied. In the present 
paper we will mostly restrict ourselves to a simple example of a fiat direction involving 
L and ("the H^L fiat direction"), which is a linear combination of superfields, and 
study its dynamical evolution through the cosmic history of the Universe. 

We begin our discussion with the basic definitions of the flat directions in = 1 
supersymmetry in Section 2. In Section 3 we first construct H^L flat direction 
parameterized by a gauge invariant monomial. We then discuss gauge invariant 
polynomials which are F and D flat at the level of unbroken supersymmetry and write 
down the equation of motion for the HuL scalars with the help of two projection 
operators, defined in the Appendix. Vanishing of the gauge currents is carefully 
accounted for. to In Section 4 we write down the multifield scalar potential, accounting 
for supersymmetry breaking and non-renormalizable superpotential corrections. In 
Section 5 we discuss qualitatively the dynamical behavior of the H^L condensate 
consisting, comparing it with single field case. Section 6 presents the results numerical 
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studies for the multifield trajectories, charge densities and charge asymmetries. In 
Section 7 we make concluding remarks and discuss the physical implications of the 
results. 



2. General features of fiat directions 



The scalar potential of supersymmetric theories is given by ^\ 



2 

A 



dW ^ 



(1) 



where W is the superpotential and the D-fiat directions are solutions to 

D^ = J2<PlT% = 0, (2) 

i 

where 0j are the scalar fields and T"^ are the generators of the symmetry group, which 
for SU{n)ca.n be chosen to be Hermitian and traceless. Making a gauge transformation 

0, ^ U{x)<p, (3) 

transforms the D-term in Eq. ^ to 

D^^Yl <PlU{xr'T^U{x)<P. = J2 <PlT% , (4) 

i i 

where T"^ is also hermitian and traceless. Since are basis for traceless hermitian 
matrices, then can be written as T"^ = c^^T^ with c^^ as complex coefficients. 
Hence the D-term, Eq. reads 

/)^^c^^D^ = 0, (5) 

where Eq. Q was used in the last step, proving that Eq. is gauge invariant. The 
solutions, Eq. (0), are gauge orbits. (This can be generalized to other gauge groups 
by virtue of Baker- Campbell-Hausdorff formula or in general quoting that D term 
transforms as a vector in the adjoint representation, which essentially have been shown 
to be the case for SU{N)). 
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One should also specify the gauge field configuration related to the scalar field 
because it is changed by a gauge transformation. The minimum energy is obtained with 
pure gauge configurations = 0. The equation of motion for the gauge field is 

D'-F^^ = , (6) 

where is the matter current defined by 

J^ = 'J2 {{D,<p,n% - 0lT^(/^,0.)) . (7) 

i 

With a pure gauge configuration the left hand side of Eq. (jH)) vanishes, providing a 
constraint on the matter current Eq. ((Zj), and a relationship between the fiat direction 
and gauge fields. Choosing a gauge corresponds to picking a scalar field value and the 
gauge field value from the gauge orbit. 

Minimizing the energy of the configurations, z.e.the equation of motion, requires 
that D^(j)i = 0. This of course results into vanishing current in Eq. ((7j). With the 
gauge choice, = 0, the scalar fields would then have to be constants. However, once 
SUSY breaking is included, the scalar fields become dynamical and the condition on the 
minimum energy becomes more subtle. For a homogeneous scalar field one can check 
that the minimum of energy is still obtained by a pure gauge configuration, despite the 
fact that the scalar fields have a non-vanishing energy. Therefore the current Eq. ((7j) 
has to vanish, but now -Do^i 7^ 0. This is a constraint on the cosmological dynamics of 
MSSM condensates that has to be satisfied. 

Let us remark that for inhomogeneous fields the situation is even more complicated. 
The minimum of energy still obeys the equation of motion Eq. (jH}, but it is no longer 
obvious that the minimum is obtained with a pure gauge configuration and a zero 
current. In fact, a more likely situation is that a spatially varying condensate, implying 
a spatially varying gauge charge, will give rise to non-vanishing currents and non- 
trivial gauge configurations. In the present paper we do not consider inhomogeneous 
condensates, but if one were to study the spatial fiuctuations of the fiat direction 
inhomogeneities would have to be taken into account. 
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The solutions of Eq. can be parameterized by analytic gauge invariant polyno- 
mials [HIIH] which can be obtained by solving 

where C is a complex coefficient |. This equation is also gauge invariant and therefore 
acts as a solution to a gauge orbit. The geometrical significance of Eq. (jH)) is that the 
polynomial / defines level surfaces, e.g. = B = const., which actually are the 

gauge orbits of the complexified gauge group. Then the flat direction corresponds to 



the minimum of the norm, |0| = 10*1^5 which turns out to be orthogonal to 0. 

For analytic functions the complex conjugate of the gradient is orthogonal to the level 
surface. This leads to Eq. (jH)) which just states that the gradient of / and the complex 
conjugate of have to be parallel (for details see 

If the polynomial describing the flat direction is actually a monomial and is 
composed of scalar fields of the model, then from Eq. (jS)) it follows that there are 
complex equations for + 1 complex variables, (pi, and C. Therefore there is at 
least one complex degree of freedom left, which can be chosen to be the scalar field 
parameterizing the fiat direction. In case of flat directions parameterized by multiple 
scalar fields, the relevant polynomial is a sum of monomials /j such that / = q/j, 
where Cj are complex coefficients §. As a consequence, there appear extra complex 
degrees of freedom. 

3. Example: HuL fiat directions 

3.1. HuL parameterized by single scalar field 

A generic treatment of polynomials generating MSSM fiat directions is beyond the 

scope of the present paper. Rather, in order to gain some insight on the dynamics of 

fiat directions with multiple fields, we focus on a simple example, the HuL fiat direction. 
I This can be generalized to the supergravity case with a non-minimal Kahler potential, K, by 
demanding dl/dc/)i — CdK/dcj)i. 

§ Actually one of the can be chosen freely, by a reparameterization of C in Eq. ||SJ . 
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Let us begin by reiterating the case for a single generation of leptons, denoted here as 
L. The relevant monomial is then given by [HI UK] 



(9) 



Applying the constraint Eq. 
dl 



dl 



, one obtains two equations 
ep^Hi = C{L^y. 



(10) 



By taking the absolute values of the equations, it follows that |C| = 1 and \H^\ 



\L^\. By multiplying the first equation by and summing over a, one 



obtains that iJ„ _L L. Hence the general structure of HuL is given by 



Hu 



\ 



(11) 



where 0, x ^-re complex scalar fields and C is a real scalar field. 

Let us now apply the current constraint = as given by Eq. ((Tj). In the case of 
single leptonic generation we can write it as 

= I {{D^HjT^H^ - HiT^D^K, + [D^LyT^L - L^T^D^L) = , (12) 

where T"^ are the SU{2) x t/(l)-generators so that there are four equations altogether. 
Since the current is gauge invariant, we can choose a gauge, and since the gauge is pure, 
we may choose = 0. (The details of the calculation is left to the next subsection. 



see Eq. 



As is well known, the final result turns out to be 

/ ,\ , /„\ 

L 



1 

71 



vV 



1 

7! 





y0y 



(13) 



where is a complex scalar field. 

One has to check that Eq. p3|) is also F-flat. The F-terms are obtained from the 
superpotential 



W = XuQHuU + XdQHdd + XeLHdC + finHuHd , 



(14) 
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where A^, A^, Ag are the Yukawa couphngs, family, colour and SU{2) indices are 
suppressed. One can easily find that the HuL is automatically F-fiat except for the 
/i//-term in the last equation. However, since hh is of the order of the SUSY breaking 
mass, and, since SUSY breaking terms anyway lift the fiat direction, the /i^-terms can 
be neglected at this stage. The Lagrangian along the H^L direction reads simply as 

C = \d,H^\^ + \d^L\^ -V= - V{ct>) . (15) 

3.2. HuL parametrized by three scalar fields 

In principle, the discussion above applies separately for each lepton generation. However, 
apart from pure chance, there is no physical mechanism that would pick out one 
generation over the others. Hence a most natural possibility would be to consider H^L 
parametrized by all three leptonic degrees of freedom. Therefore we should consider the 
gauge invariant polynomial 

/ = viHuLi + U2HuL2 + usHuLs , (16) 

where z/j are complex coefiicients (of which one can be freely chosen). The D-fiatness 
equations Eq. (jSJ become (compare with Eq. (UHl)) 
dl 



OH" 
dl 



i=l 

i^.e^^H^ = C{Lty . (17) 



dLf 

By solving from the second equation and inserting into the first, one obtains the 
constraint 



|2 



Working the other way around, i. e. by solving iJ" from the first equation and inserting 
into the second, one obtains a matrix equation 
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where L° = (L", , L'^y" is a vector in flavor space and P is a projection matrix to the 
vector V = [ui, 1/2, 1^3)^- Note that Eq. |T!Hl can be fulfilled only if (L")* is parallel to u 
so that 



(20) 



where is a complex coefficient. The fiat direction is given by 

= V. 



„2* 



(21) 



where Sc is the phase of C which was left undetermined in Eq. fll8j) . Making redefinitions 
(pi = z/*c^*, ip = c^*/c^* and 6c = X ~ ^2, where x is real field and 62 is the phase of c^, 
one obtains a simplified expression for the field content of the if^L-fiat direction 

\ 



-r 



(22) 



Let us now apply the current constraint = Eq. ((7)), where 
= i{{D,H^)^T^Hu - HlT^D.H^) + z 5^((D^L,)tT^L, - llT^D.Li] 

i 

We find 



J, 



where 



^ i fc / k 

Y^jf-2tY^\cp,\^d,x \ - 2 5^ 10,1 

^ i fc / k 



Adding Eqs. ^IM 1^ together, we obtain 



(23) 

(24) 
(25) 
(26) 

(27) 
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Using Eqs. (^Sl YI^ . one obtains d^ip = 0, so that ip = const. Therefore we may gauge 
it away by a global U{1) transformation including it in the normalization of (pi. Hence 
we finally obtain 



_ (29) 

where x is given in Eq. (j28j) . Since the final configuration is three complex dimensional 
which is the maximal dimension of the fiat direction J3] containing and Li, the 
polynomial Eq. gives a complete characterization of the HuL fiat direction. 

3.3. The effective Lagrangian for the multiple flat direction 

For one scalar field the constraint for the phase x in Eq. (j2HI) reduces to d^x = 9^9, 
where 6 is the phase of the scalar field = |0| exp(i^). Therefore, up to a normalization, 
the field configuration reduces to Eq. !^l'd\i. However, for more than one scalar field, the 
constraint Eq. ()28j) is not necessarily solvable. The problem is the following: on the 
right hand side there is a vector field (formed as a sum of the scalar fields (pi and their 
derivatives), on the left hand side there is a gradient of the scalar function, so the right 
hand side has to be rotationless (in the language of differential forms the one-form on 
the right-hand side has to be closed for it to be exact). The necessary (and in simply 
connected space sufficient) condition for the vector field to be a gradient of a scalar 
function is 

{d^d,-d,d^)x = 0, (30) 

which yields 

9^$t (1 - Pi -P2) 5.^ = 0, (31) 

where the Pi are projection operators discussed in the Appendix. 

From Eq. ()31|) one sees that the components of the gradient are constrained to 
be orthogonal outside of the two-dimensional surface spanned by $ and \E'. On the 
surface there are no constraints. If there is only one scalar field, then $ and \E' are 
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two-dimensional and Eq. (j31|) is trivially satisfied, since in this case Pi + P2 = 1. If the 
scalar fields do not depend on space or time, then the constraint is trivially satisfied. 
If the fields depend only on one coordinate (time or one of the spatial coordinates; 
homogeneous fields fall into this category), then again one can show that the constraint 
is trivially satisfied. Since in the cosmological context the large scale homogeneity is a 
reasonable approximation, we do not have to worry here about not fulfilling Eq. ()31|) . 

The dynamical evolution of the multiple field fiat direction is determined by the 
MSSM Lagrangian 

3 

C=\d^Hu\'' + J2\W-y^ (32) 



i=l 



It turns out to be convenient to change the parameterization of the flat direction by a 
phase shift to 



where x is solved from Eq. 
into Eq 




(33) 



Now the Lagrangian is obtained by inserting Eq. ()33p 



C = -9^$^ (^1 + Pi - -P2 j df'^-V, (34) 

(for the projection operators Pi, see Appendix). Note that there are now non- minimal 
kinetic terms. This is due to the fact that the HuL direction is a curved sub-manifold 
of the whole field manifold, as can be seen from Eq. (jH^ which implies that, up 
to a gauge choice, the sub-manifold is actually a sphere. In the one-field case the 
fiat direction is formed only along a one-dimensional sub-manifold and therefore has 
vanishing curvature. The equations of motion resulting from Eq. (jH^ are (see Appendix 
for details) 



dV 



+ ^ {df,<I/^P2d''^) + -$ 9^.$^ ( 1 - Pi - -P2 I df"^ 



(35) 



where = 
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Since we are interested in the background dynamics, the partial derivatives can be 
replaced by time derivatives. In that case the different terms have analogues in classical 
mechanics. On the first row the potential gradient is deformed by two projections: Pi 
makes the potential flatter and P2 steeper in the directions $ and \1/ respectively. The 
last term on the second row generalizes centripetal acceleration. The rest two terms are 
analogous to the Coriolis and Euler forces. 



4. The potential for HuL 



In the real world supersymmetry is broken. This lifts the degeneracy of the vacuum 
solutions of the supersymmetric gauge theory, and the fiat directions become dynamical 
fields 13 in]. Here we just list the contributions relevant for H^L. 

The soft SUSY breaking mass parameters together with the //-term read 

{rnl^+t^li)\Hu\' + J2<\L^\" = Y.<\^^\"^ (36) 

i i 

where 



2 2,2,2 



(37) 



The mass parameters m^^ are of the order of the soft SUSY breaking mass nis, which 
in the gravity mediated SUSY breaking scenario is the gravitino mass ^3/2- The non- 
renormalizable operators lifting H^L are given by fl^ . 

^ ^ij i^uLi) (H^Lj) , (38) 



W 



where Xij = Xji are (complex) coupling constants and M is a large mass scale (typically 
GUT or Planck scale). The superpotential Eq. (|38p gives rise to potential terms through 
the F-terms 



dW 

1 

4M2 



i 



dW 



(39) 
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as well as through the generalized A-terms 

Am^/^W + h.c. = I'^'^l' H ^^^'^^'^^ + ^-^^ ' (40) 

k ij 

where A ~ 1 is a complex coefficient. If the Kahler potential contains a non-trivial 
coupling to the inflaton, such as 

K~-^/t/0V^ (41) 
where / is the inflaton and is if„ or Lj, then there is a mass correction 

Y.c,H^\H\ (42) 

i 

where H is the Hubble parameter and |cj| ~ 1 after inflation. During F-term inflation 
|cj| ~ 1, but for D-term inflation |cj| ^ 1 20j. Usually we take Cj = —1. There are also 
Hubble induced A-terms if the inflation is induced by the F-term 

aEW + h.c. , (43) 

where a ~ 1 is a complex coefficient, for D-term inflation |a| -C 1 |2nj . 

5. Dynamics of the multifield H^L fiat direction 

5.1. Motion of HuL during inflation 

During inflation H ^ Hj = const, and the minima of the potential are fixed points of 
the equation of motion for the homogeneous mode Eq. (|35|) as in This follows from 
the fact that $ = $ = is a fixed point so that eventually only the potential term in 
Eq. ()35|1 remains in the equation of motion. Local stability is due to the fact that the 
last three terms of Eq. p5|) do not contribute to linearized perturbations. 

5.2. Motion of HuL after inflation 

After inflation the Universe is dominated by the inflaton oscillations, which produces 
effectively a matter dominated Universe with H = 2/(3t). Before the oscillations start, 
H ^ nis, we can make the approximation that only the Hubble induced mass 
term Eq. (j42j) and the non-renormalizable terms Eq. (|39p are important. (The Hubble 
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induced A-term Eq. (|43p affects only the phase). Then we can find a fixed point solution 
following [Hj and making the change of variables 

(t>, = {MHfl^i,,, t = m^'e'. (44) 

At the fixed point, dxi/dz = 0, so that the equation of motion simplifies to 

where Vnr is the potential in Eq. with ipi as independent variables instead of 
and C is a matrix with Cij = Ci6ij. An effective potential can be obtained for Eq. 
whose extrema are fixed points of the equation of motion. This is a tedious calculation 
in general but with q = c for all i one finds 

K// = M'V^nii;) -2(^ + 1] • (46) 



. 9 4, 

Hence the effective potential is bounded from below provided Vatj? is. The minima of 
Veff correspond to marginal fixed points of the equation of motion. This reasoning is 
similar to the one field case ,6j, where there is no friction term for d = 4 fiat directions. 
However, one cannot solve Eq. in a closed form, but solution(s) definitely exist since 
Vnr is positive definite and higher order than the second term. When all the terms in 
the potential are taken into account, this is only an approximate behavior. 

5.3. The onset of rotation H ~ nis 

hi H ms the scalar fields </)j feel the torque produced by the A terms and start to 
rotate. The precise onset of rotation depends on m^. and Cj because the phase rotation 
of (pi starts when \ci\H'^ ~ m^.. When H <^ nis the lepton charge densities in the 
co-moving volume of 0j asymptote to (different) constant values. The scalar fields at 
this time decay as a{t)^^/'^ for a matter dominated Universe. This behavior is similar to 
the one-field case IH] • However, for quantitative results one has to resort to numerics. 
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Figure 1. The trajectories of (pi for a — i, Ci — A = —1 and m^. — lOOGeV with 
dashed line for the first, dash-dotted for the second and dotted for the third family. 



6. Numerical results 



Since we expect a fixed point behavior, for numerical calculations it is useful to choose 
the variables as defined in Eq. ()44|) . In the one-field case it is easy to obtain the initial 
conditions for the evolution after inflation by solving for the field value at the minimum. 
Here this is not possible. Instead we assign random initial values and solve the equations 
of motion during inflation numerically for 50 e-foldings and then insert the values thus 
obtained into the equations of motion valid after inflation. In practice the duration of 
inflation is inessential as the fields tend to relax into the minima in 10-20 e-folds. 

For simplicity, we have restricted our considerations to the case where the Yukawa 
coupling matrix Xij is diagonal. We then find inumerical solutions which are stable. In 
Figs. H 121 we have plotted the trajectories of the scalar fields (pi in the complex plane. 
In contrast with the one field case [01121! the scalar fields do not relax directly into the 
origin when H <^ms. This is also implicit in Fig. El where there is a non-zero charge 
for the slepton fields right from the beginning. This feature can be ascribed to a phase 
oscillation around the minimum of the potential. 

In Fig. 121 we have plotted the evolution of the lepton charge densities in the co- 
moving volume of the slepton fields Li as parameterized in Eq. ()33|). Note that the total 
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Figure 2. The end trajectories of (pi for a = i, Ci = A = ~1 and m^. = lOOGeV with 
dashed hnc for the first, dash-dotted for the second and dotted for the third family. 



lepton charge is in units of m|Mp, where = 100 GeV. The initial conditions were 
again chosen by first solving the equations of motion during inflation. 

In Fig. m we have plotted the dependence of charge density against the phase 
6a of Hubble-induced A-term coefficient a with two different initial conditions during 
inflation. The different initial conditions definitely affect the charge densities of the 
individual slepton fields. However, although the initial conditions also affect the total 
charge density, the effect is not as large as for the individual charges. If one where to 
choose all the initial field values equal so that all the charge densities would be equal, the 
total charge density evolution would differ from those presented in Fig. |3]in magnitude 
but not in overall shape. 

In Fig. El we plot the charge asymmetry qtot/ntot, where qtot = and Utot = 
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family. The total slepton charge density is given by the solid line. The parameters 
were m^. = 1, Ci = A = —1, a = i and Xij = 6ij. 
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Figure 4. The charge density vs. the phase Oa of a with initial conditions during 
inflation = 1, ^2 = i and (f)^ = —1 in the upper figure and = 1, ^2 = « and 
(j)3 = —1 + i in the lower figure, where field values are in units of {MpHiy/^. The 
initial velocities vanish and the parameter values are m^. =1, Cj=^ = — 1, \a\ = 1 
and Xij = 5ij. Again solid line is for total charge, dashed line for Li, dash-dot for L2 
and dotted for L3. In the upper figure Li and L3 are actually the same curve. 
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Figure 5. The charge asymmetries related to the cases of fig. O with the same 
parameter values and initial conditions. 



^ \qi\. One sees that for all these cases one produces a charge asymmetry of roughly 
~ 0.1. 

7. Conclusions and Discussion 

The numerical studies indicate the produced baryon-to-entropy ratio is not significantly 
changed when including all the lepton families to the HuL fiat direction. It thus 
seems that the simple one-field case captures the essential features of the Affieck-Dine 
mechanism, at least for H^L. 

However, there are some interesting details that do change. The first is that the 
flat direction has non-minimal kinetic terms, which is due to the fact that the sub- 
manifold spanned by the flat direction is curved. In the one-fleld case non-minimality 
is absent because all the one- dimensional manifolds are flat. Hence there is always a 
choice of coordinates where the fleld metric is given by just the Kronecker delta. Another 
difference is that there is phase motion of the flelds already right after inflation. This 
is due to the F-terms containing cross terms that break individual lepton symmetries. 
However, the F-terms conserve the total lepton charge which is produced only through 
the A-term. The individual charges may in general have opposite signs and can give 
substantial contribution to the entropy, which is proportional to the sum of absolute 
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values of the individual charge densities. The total charge is the sum of individual 
charge densities and thus cancellations can occur. 

Although the behaviour of a homogenous H^L background with multiple scalar 
fields is not qualitatively different from the corresponding one-field case, this may not be 
true if inhomogeneities are taken into account. For small cosmo logical perturbations and 
multiple fields the gauge current condition of Eq. (jHT|) becomes highly non-trivial. This 
can be seen as follows. HuL is described by three complex scalar fields or six real degrees 
of freedom. The current condition Eq. ()31|) defines a two-dimensional submanifold in 
the six-dimensional space, so that there are four real fluctuation degrees of freedom that 
would not obey Eq. (jHlll . Hence there would exist matter currents. This implies that 
the pure gauge field no longer provides a solution for the equations of motion. As a 
consequence, there would arise electric and magnetic gauge fields which could contribute 
to the CMB perturbations. These might be of paramount importance for condensate 
fragmentation, since the fragments (the ground states of which are Q-balls) result from 
the non-linearities of the perturbations [3 |Hj. The Q-ball properties PX] have been 
considered only for condensates with minimal kinetic terms. 

In this paper we focused on the H^L direction. For other directions more 
complications may arise. For example, the d = 6 fiat direction {LL){ddd){udd) has six 
different superpotential terms that contribute to the lifting of the fiat direction: d = 4 
term uude, d = 5 terms HuLudd and HdLddd and d = 6 terms uddQLd, {udd){udd) 
and LLeudd. The c? = 5, 6 terms have B — L = —2 whereas the d = A term has 
B — L = 0. These will give rise to cross terms in the F term which has a non-zero 
B — L and therefore contributes to the total charge. Such a situation occurs whenever 
the various superpotential terms with different B — L have at least one field in common. 
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Appendix A. Projection Matrices 



Projection of a vector a along a vector b is given by 

, < bla > / A 1 N 

b^^, (A.l) 

where < -I- > is a (complex) scalar product. The operation of projection corresponds 
to a matrix 

bbt 



P 



btb 



(A.2) 



where we have marked the vector with a column matrix and f represents the hermitian 
conjugation. Then Eq. ()A.1|) is just Pa. 

In the present paper the projection operators needed are somewhat complicated 
since in principle they operate on the real basis of complex vectors. Therefore we write 
complex scalar fields in terms of two real fields 0j = ipi + iXi- Then we can form a vector 



(A.3) 



which represents in the real basis. One can change the basis to complex vectors and 
(h* with 



/ 











V 


r) 









and construct a projection operator 
$$t 



(A.4) 



Pi 



We also need a projection operator along the vector orthogonal to $ 



P9 



where 



v[/tvl/ ' 



(A.5) 



(A.6) 



(A.7) 
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It is easy to check that = p2 ^ ^^^^ p^p^ ^ p^p^ ^ q ^j^j^ ^.^g ^gjp ^.j^gg^ 

projection operators we can form the following combinations 



1 2 



2E 



1 2 



2E 



With these identities we can write the Lagrangian as 



£ = ^5^$^ (^1 + Pi - ip2^ - Vi<^, <l>t) . 



(A.8) 



(A.9) 



The equations of motion read then 



1 \ ■ dV ^t^M^ 

1 + A - (9„a-* + 3i/<i.) + - a,*^ 



-^a„<i>t ( 1 - Pi - -P2 1 = , 

<i>t<i> ^ V 2 7 



(A.IO) 



$t$ A* " " ^"^ " ' <l,t<l,^ 
where a is the scale factor of the FRW metric. 

There is a formula for inverting a unit operator plus a sum of orthogonal projection 
operators, given by 



-p. 



(A.ll) 



where P^P,- = 5ijPi and G C\{-1}. Using Eq. ()A.11|1 in Eq. ()A.10|1 . we then obtain a 
form for the equation of motion which is more easily comparable to the standard case: 



^ 1 „ ^\ dV ^ ^ 
^-2^^ + ^0a$t-^^^W 



^d,¥P,d^<^ + ^a,vl.t (1 _ _ 3p^^ = . (A.12) 
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